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The problem of precise evaluation of perturbative QCD predictions at moderate energies is addressed. In order 
to improve stability of the predictions with respect to change of the renormalization scheme it is proposed to 
replace the sequence of conventional renormalization group improved approximants by a sequence of modified 
approximants, involving a modified running coupling constant, which is free from Landau singularity and much 
less renormalization scheme depenedent than the conventional running coupling constant. A concrete model of 
the modified coupling constant is proposed and it is shown, that the QCD corrections to the static interquark 
potential evaluated with this coupling constant are indeed much less sensitive to the scheme parameters. It 
is pointed out, that the modified predictions display somewhat weaker energy dependence compared to the 
conventional predictions, which may help accomodate in a consistent way some very low and very high energy 
determinations of the strong coupling constant. 


As is well known, finite order QCD predictions 
obtained with the conventional renormalization 
group improved perturbation expansion depend 
to some extent on the choice of the renormaliza¬ 
tion scheme (RS) pp. This scheme dependence is 
formally of higher order relative to the order of 
the perturbative approximant, but at moderate 
energies (of the order of few GeV) it may be nu¬ 
merically significant |5]. This casts some doubt 
on the accuracy of perturbative results obtained 
in this energy range. Some attempts have been 
made to address this problem by modifying con¬ 
ventional perturbative approximants [^. In the 
following we present an alternative approach 
which may offer some advantages over the previ¬ 
ously proposed solutions. 

The starting point of our analysis is the ob¬ 
servation, that strong RS dependence of pertur¬ 
bative predictions is to a large extent the result 
of a very strong RS dependence of the running 
coupling parameter (couplant) itself. The N-th 
order couplant a(Q^) = g^(Q^)/47r^ in the con¬ 
ventional approach satisfies the renormalization 
group (RG) equation: 

( 1 ) 


with 
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where the coefficients C2, C3... are scheme depen¬ 
dent. An extreme manifestation of the strong RS- 
dependence of a{Q'^) is the fact, that in a large 
class of renormalization schemes (in particular, 
all the schemes with positive Ck coefhcients) the 
couplant becomes singular at positive (Lan¬ 
dau singularity), and the location and charac¬ 
ter of the singularity strongly dependends on the 
scheme parameters (it is controlled mainly by the 
highest order term retained in the /3-function). 

However, within the perturbative approach we 
have some freedom in defining the effective expan¬ 
sion parameter. The idea is to exploit this free¬ 
dom and define an alternative, less RS-dependent 
couplant, which could then be used in perturba¬ 
tive expressions instead of the conventional cou¬ 
plant, hopefully giving less RS-dependent predic¬ 
tions for physical quantities. Obviously, we would 
like this modified couplant to be free from the 
Landau singularity. 

A simple way to define such a modified cou¬ 
plant is to replace the conventional, polynomial 
/3-function @ by an appropriately chosen non- 
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polynomial function /3(a). To retain the formal 
consistency with the conventional approach, the 
expansion of the modified /3-function in any given 
order should of course reproduce the required 
number of small-a expansion coefficients of the 
conventional /3-function. Next, to ensure the ab¬ 
sence of the Landau singularity it should behave 
for large a as —^a^, with ^ > 0 and fc < 1. Fur¬ 
thermore, it would be convenient, if the modified 
/3-function would be analytic in some neighbour¬ 
hood of a = 0. This restriction eliminates mod¬ 
els with exponentially small contributions, which 
produce 1 /Q” terms in large-Q^ expansion for the 
couplant and the physical quantity itself, that in¬ 
terfere with terms of similar type from the oper¬ 
ator product expansion, and sometimes are even 
forbidden in this expansion. 

Additionally, we shall impose the condition, 
that (a) should not contain any order-specific 
free parameters, i.e. it should contain only pa¬ 
rameters, which characterize the whole sequence 
of the modified approximants. Indeed, if (a) 
contains in each order new free parameters, un¬ 
related to other parameters, this would have ad¬ 
verse effect on the predictive power of the modi¬ 
fied expansion, at least if we restrict our attention 
to few low order approximants. 

The modified couplant is obtained by integrat¬ 
ing the renormalization group equation o with 
the modified /3-function. 

The idea of modifying the expression for the ef¬ 
fective coupling parameter in QCD has of course 
a long history, and many models have been pro¬ 
posed in the literature 0. It turns out, however, 
that these models are not satisfactory from the 
point of view described above. 

A simple way to construct models of 
with the required properties is to use the so called 
mapping method . One of the models analyzed 
by us is based on the mapping 


where 77 is a real positive parameter. It is easy to 
verify, that the function of the form 
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k=0 

where At is a real positive parameter and the co¬ 
efficients Ck have the form 

Co = 1 — r]K, Cl = Ci + 2r] — 

C2 = C2 + 3ci77 -I- 3?7^ — ry^K, etc. 

does indeed satisfy all the criteria listed above. 

The free parameters in this model function 
have the following interpretation: \/rj character¬ 
izes the value of the couplant a(Q^), for which the 
nonpolynomial character of the / 3 *^^^ (a) becomes 
essential, while k determines the low-Q^ asymp¬ 
totics of the modified couplant a{Q^), which be¬ 
haves in this limit as The func¬ 
tion is of course a viable replacement for 

/ 3 W(a) in the N-th order of perturbation expan¬ 
sion for any value of these parameters. However, 
a clever choice of these parameters, based per¬ 
haps on some information outside of perturbation 
theory, may improve the quality of low order per¬ 
turbative results in the modified expansion. As a 
first try, we choose k = 2/6, which ensures 1/Q^ 
behavior for the modified couplant at low Q^, as 
suggested by some theoretical approaches. Our 
procedure for fixing i] is described further below. 

The modified perturbative predictions for a 
physical quantity are obtained by replacing the 
conventional couplant in the perturbative expres¬ 
sion by the modified couplant. To verify, whether 
our modified couplant does indeed improve the 
stability of perturbative predictions, we analysed 
in some detail the next-to-leading (NL) and next- 
next-to-leading (NNL) order predictions for sev¬ 
eral quantities, including the QCD correction dy 
to the static interquark interaction potential 0, 
which is related to the Fourier transform of this 
potential in the following way: 

vm = ( 5 ) 

where Cp =4/3. The NNL order expression for 
i5v has the form: 

^ (Q^) = « [1 + ri a -I- r2 a^] , 


( 6 ) 
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For three active flavors in the MS scheme with 
^2 = q 2 have = 1.75 and = 16.7998 
(as well as 5 = 9/2, a = 16/9 and c^® = 4.471). 
The freedom of choice of the RS for this quantity 
in NL order may be parametrized by the parame¬ 
ter ri, while in the NNL order we have additional 
free parameter, which we may choose to be the 
coefficient C 2 in the /3-function. The numerical 
value of the couplant is then determined from the 
implicit equation, which results from the integra¬ 
tion of the of the equation m with an appropriate 
boundary condition: 
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In Fig. n we show predictions for Sy at = 
3GeV^, as a function of ri, for several values of 
C 2 , obtained in the conventional expansion. As 
we see, the differences in the predictions are quite 
large. In Fig. [21 we show the corresponding plot 
obtained in the modified expansion (with r] = 4.1, 
which we justify further below). Clearly, the 
modified predictions are much less RS-dependent 
than the conventional predictions. We have 
verified, that similar improvement in the stabil¬ 
ity is obtained for other physical quantities, in 
particular for the QCD correction to the Gross- 
Llewellyn-Smith sum rule in deep inelastic scat¬ 
tering |5|. 

The pattern of ri and C 2 depedence of the mod¬ 
ified predictions indicates, that also in the case of 
the modified expansion it is possible to choose 
the scheme according to the Principle of Minimal 
Sensititivity |2]. General equations defining the 
PMS parameters for the modified expansion are 
of course the same as in the case of the conven¬ 
tional expansion (vanishing of the partial deriva¬ 
tives with respect to the scheme parameters), but 
the resulting equations are much more compli¬ 
cated, because the integral of l//3(a) is more in¬ 
volved. Let us note, that it is only in the modi¬ 
fied expansion that the full potential of the PMS 
method may be realized, because only in the mod¬ 
ified approach the predictions are finite down to 
= 0 in all possible renormalization schemes. 
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Figure 1. 5y for uf = 3, at = 3GeV^, as a 
function of ri, for several values of C 2 , as given 
by the conventional perturbation expansion with 
= 350 MeV. Dashed line indicates the NL 
or^er prediction. 


An important point in our analysis is the choice 
of the value of the parameter ry. To fix the value 
of r] we use the phenomenological formula cni 
for i5v, which has had some success in correlating 
the experimental data for heavy quarkonia. We 
adjust ry in such a way, that the dependence 
of the NNL order modified PMS prediction for 
6v would match the Q^-dependence of the phe¬ 
nomenological expression m as closely as pos¬ 
sible. It turns out, that the best agreement is 
obtained for 7y = 4.1 — if the modified predic¬ 
tions for Sv are matched to coincide with the phe¬ 
nomenological expression at = 9GeV^, then 
the relative difference between these two expres¬ 
sions is less than 1% down to = 1 GeV^. 

Perturbative predictions obtained in the modi¬ 
fied approach have an interesting property: they 
have somewhat weaker Q^-dependence at moder¬ 
ate compared to the conventional predictions. 
In order to see, how significant this effect could 
be for phenomenology, we performed within the 
modified approach a fit to the experimental data 
for ^gls- From the experimental result given in 
m we inferred, that purely perturbative contri¬ 
bution to the GLS integral at 3 GeV^ is equal to 
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Figure 2. Sy for uf = 3, at = 3GeV^, as a 

function of ri, for several values of C 2 , obtained 

using the modified perturbation expansion with 

A— = 350 MeV and n = 4.1. Dashed line indi- 
MS ' 

cates the NL order prediction. 


= 0.131 ± 0.040, where we have added the 
statistical and systematic errors in quadrature. If 

/q\ 

we fit to this experimental value using the 
conventional NNL order expression for i5gls in 
the MS scheme, and then convert this parame¬ 
ter into the corresponding value of us¬ 

ing the matching formula m to cross the quark 
thresholds, we obtain q;s(M|) = 0.1126tooii6i 
which is to be compared with the the world av¬ 
erage as(M|) = 0.1187±0.0020 Performing 
the same fit and extrapolation with the modified 
NNL order expression in the MS scheme (with 
r] = 4.1) we obtain a somewhat larger value: 
as{Ml) = 0.11471o;q 526- ^he upward shift is a 
welcome effect. This effect is even more strik¬ 
ing, when we perform the same fit using the PMS 
approximants. Using the NNL order PMS ex¬ 
pression in the conventional expansion we obtain 
q:s(A 7|) = 0.1097to'oio2i while using the NNL or¬ 
der PMS approximant in the modified expansion 
we get ds(A7|) = O.llSOto^oi*?- This shows, that 
within the the modified perturbative approach it 
may be easier to accomodate in a consistent way 
some of the very low and very high energy deter¬ 
minations of as. 
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